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Particle production of an Abelian vector boson field with an a:xial coupling is investigated. The 

conditions for the generation of scale invariant spectra for the vector field transverse components are 

obtained. If the vector field contributes to the curvature perturbation in the Universe, scale- invariant 

fN^ ■ particle production enables it to give rise to statistical anisotropy in the spectrum and bispectrum 

of cosmological perturbations. The axial coupling allows particle production to be parity violating, 

C^ , which in turn can generate parity violating signatures in the bispectrum. The conditions for parity 

^Sl ■ violation are derived and the observational signatures are obtained in the context of the vector 

_^ ' curvaton paradigm. Two concrete examples are presented based on realistic particle theory. 
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I. INTRODUCTION 

In the last few years the contribution of cosmic vector fields to the curvature perturbation in the Universe is under 
investigation. This effort was triggered by the pioneering work in Ref. [l|, which introduced the vector curvaton 
mechanism, through which vector boson fields can contribute to or even fully generate the curvature perturbation 
(for a recent review see Ref. Qj). In analogy with the scalar curvaton mechanism [3|, the vector curvaton is a 
spectator field during cosmic inflation, which becomes heavy after the end of inflation, when it can dominate (or 
nearly dominate) the density of the Universe before its decay, thereby imprinting its contribution to the curvature 
perturbation C,. It was soon realised that, through their contribution to (, vector flelds can give rise to distinct 
observational signatures, namely produce statistical anisotropy in the spectrum and bispectrum of ( [^-[Sl- This is 
a new observable, which amounts to direction dependent patterns on the CMB temperature anisotropics, which can 
be predominantly anisotropic in the bispectrum [a, 01 . Such patterns cannot be generated if scalar fields alone are 
responsible for the formation of (^. Thus, in view of the imminent observations of the Planck satellite mission that 
may well observe statistical anisotropy, it is of paramount importance to investigate this effect and what it can reveal 
to us for the underlying theory. In particular, if statistical anisotropy is observed by Planck, we will have our first 
glimpse into the gauge field content of theories beyond the standard model. 

In general, one can divide the mechanisms that have been developed so far for the contribution of vector fields into 
^ in two classes. Firstly, one can have an indirect infiuence of the vector field to inflation. One way to do this is by 
considering that the vector field develops a condensate, which generates anisotropic stress that renders inflation mildly 
anisotropic. In this case, the anisotropy in the expansion reflects itself onto the perturbations of the inflaton field, 
thereby generating statistical anisotropy in C [1?;TJ| (for a recent review see fl4!|). This anisotropisation of inflation 
is usually achieved by introducing some coupling between the vector field and the infiaton [12 - 15] . Alternatively, such 
coupling can backreact to the generation of the inflaton perturbations regardless of anisotropising the expansion [1^1 ■ 

The other class of mechanisms which give rise to statistical anisotropy considers the direct contribution of vector 
field per turbations to C. One way this can be done is through the vector curvaton mechanism mentioned above 
Uli M, [17h21| . where perturbations of the vector field perturb its energy density and, hence, the moment of (near) 
domination of the Universe by the vector field condensate. Statistical anisotropy in this case is due to the fact that 
vector fields undergo anisotropic particle production, in general [a, [l^. Other mechanisms have been employed as 
well, such as the end of inflation mechanism [J, |5|. Finally, one way to use directly the perturbations of vector flelds 
to source C is by considering a large number of them acting as inflatons as in Ref. |22i| . 

What do observations say about statistical anisotropy in the curvature perturbation? The observed bound on 
statistical anisotropy in the power spectrum of ( is surprisingly low; it seems that as much as 30% of it is still allowed. 
In fact, statistical anisotropy at this level was reported in Ref. [10] at the level of 9-(jI However, the direction of 
the anisotropy was suspiciously close to the ecliptic plane so the authors of these studies conclude that the flnding is 
probably due to some systematic mistake, hence being treated as an upper bound. The Planck mission will reduce this 
bound down to 2% if statistical anisotropy is not observed [9]. This means that the so-called anisotropy parameter, 
which quantifles statistical anisotropy in the spectrum of C, must lie in the following range if it is to be observed in 
the near future: 

0-02 < 5(; < 0.3 . (1) 

Although, in principle, statistical anisotropy in C can be scale dependent, a scale invariant spectrum of vector 
field perturbations results in scale- independent (7^, which does not have to be fine-tuned such that it falls into the 
above range on the cosmological scales. Moreover, even if Eq. ([1]) is satisfied for cosmological scales, a strongly tilted 
spectrum would generate intense anisotropy, giving rise to excessive curvature perturbations (e.g. leading to copious 
primordial black hole formation) or destabilising inflation itself. Thus, a (nearly) scale-invariant spectrum of vector 
field perturbations is preferred. As shown in Ref. [13J , this can be naturally attained in certain types of theories. 

Statistical anisotropy generates angular modulation of the power spectrum, but also of higher order correlators. 
Moreover, even if the anisotropy in the spectrum satisfies the upper bound in Eq. ([T]), higher order correlators can be 
predominantly anisotropic (see e.g. [19]). It is thus important to study the effects of statistical anisotropy on hi ghe r 
order correlators ^gj] as well as develop methods of detecting these effects in the CMB temperature perturbation |Tl| . 

Not many models exist as yet, for the formation of a super horizon spectrum of vector field perturbations during 
inflation. The problem is rather old as it was investigated, at first, in order to generate a primordial magnetic field 
during infiation, with superhorizon coherence. A massless Abelian vector boson field, such as the photon, cannot 
undergo particle production during inflation because it is conformally invariant. A breakdown of its conformality is, 
therefore, required in order to form the desired superhorizon spectrum of perturbations. Modiflcations of the theory in 
order to attain such breakdown were originally investigated in Ref. [23], by coupling electromagnetism non-minimally 
to gravity. Some of these proposals were recently implemented in the effort to generate a contribution of vector field 
to ^; notably a non-minimal coupling to gravity of the form RAf^A^. Such coupling was employed in Ref. J22l| where 



hundreds of vector fields are used as inflations, and also in Refs. p, [l^l, in the context of the vector curvaton model. 
However, this proposal was criticised for giving rise to ghosts |2J| (see however Ref. |l25]). 

Another model, which does not suffer from instabilities, is considering the supergravity-inspired varying kinetic 
function for the vector field fltjFuj^F^'^ ■ The latter has been shown to give rise to a new inflationary attractor 
under fairly general conditions |12l.ll3|. which leads to scale-invariant vector field perturbations [13J , when the kinetic 
function is modulated by the inflaton field. This model also has a long history, since it has been used to generate 
primordial magnetic fields too [26]. More recently, however, it has been used to affect ( in the vector curvaton 
mechanism [5, 17, 19-21)] as well as the end of inflation mechanism |4l.l27fl. 

A third popular choice for primordial magnetic field generation during inflation is considering the axial term 
h{t)F^i,F^^'^ [281. Recently, some works have investigated the vector field backreaction onto inflation, if the axial 
coupling h is modulated by a pseudo-scalar inflaton field. Ref. [20| suggests that the axial term, can allow steep 
infiation with sub-Planckian axion decay constant, evading thereby the basic problem of natural inflation [30]. In 
Ref. Ll^], the backreaction of the generated vector fleld perturbations onto the inflaton perturbations is shown to 
generate signiflcant non-Gaussianity in the latter. Finally, the effects of the axial coupling onto gravitational waves 
is investigated in Ref. [31], which may provide parity violating signatures on the tensor modes and their effects onto 
the CMBQ 

Parity violation is a special property of the axial model, compared to the other two mentioned above {fF'^ and 
RA^), which are parity conserving. What would it amount of if the vector field perturbations directly contributed 
to ^, e.g. through the vector curvaton mechanism? As shown in Ref. [5]] parity violation does not feature in the 
power spectrum of the density perturbations. However, in Ref. [6| it was found that parity violation does affect the 
bispectrum of ^ and it should reflect itself onto the shape and amplitude of the anisotropic /nl • In this paper we 
present a complete study of the particle production of an Abelian vector boson fleld featuring both a varying kinetic 
function and a non-trivial axial couplingjj We obtain scale invariant parity violating spectra and we single out the 
conditions for their successful generation, providing some concrete examples based on particle theory. Finally, we 
apply our findings onto the vector curvaton paradigm and find the generated g(^ and /nl, which are to be contrasted 
with observations. 

The structure of our paper is as follows. In Sec. [H] we present a brief review of the general case when a vec- 
tor field contributes directly into the spectrum and bispectrum of ^. In Sec. IIIII we present the axial model and 
investigate particle production throughout the parameter space. In Sec. IIVI we focus on the most promising case 
which produces scale-invariant parity violating spectra for the vector fleld components and apply our findings to 
the vector curvaton mechanism. In Sec. |V] we discuss two concrete examples based on realistic particle theory. We 
conclude in Sec. IVIl Throughout our paper we consider natural units where c = h = ks = ^ and Newton's grav- 
itational constant is SttG — rrip^, with mp = 2.4 x 10^® GeV being the reduced Planck mass. We use the metric 
Qfii^ = diag (l, — a^, — a^, — a^) and assume (quasi) de Sitter inflation with Hubble parameter H « constant. 

II. THE ANISOTROPIC, PARITY VIOLATING /nl 

In this section we briefly summarize the results in Refs. [5|, |6| and calculate the anisotropic /nl in the so called 
"flattened" shape for the first time. 

Massive vector fields have three degrees of freedom. During inflation, after horizon exit quantum fluctuations 
of these degrees of freedom become classical perturbations. To deal with these perturbations it is convenient to 
use circular polarization vectors. As these vectors transform differently under rotations, Lorentz invariance of the 
Lagrangian guarantees that equations for the three polarizations are uncoupled. Let us denote the power spectrum of 
left- and right-handed polarization modes by Vh and Vb, respectively and the longitudinal one by V\\ . Furthermore, 
expressions for the power spectrum and higher order correlators of the primordial curvature perturbation ^ become 
much simpler expressed in terms of V± , which are defined as 

V± = ^ {Vk ± Vi.) , (2) 

where V+ correspond to parity conserving and V- parity violating spectrum. The later is non-zero only if the left- 
and right- handed polarizations acquire different perturbation spectrum. Calculating /nl we will also find convenient 
to normalize spectra as 

^ other works on parity violation in the graviton bispectrum can be found in Ref . |32||. 
^ The contribution of non- Abelian vector fields to f has been considered in Refs. |27Ll34l|. 



If the p and/or q parameters are non-zero the vector field perturbation is statistically anisotropic. Note that by 
definition the p parameter can take values p> —\ and the parity violation one — 1 < g < f . 

To calculate the curvature perturbation C, we use the so called 5N formalism. This formalism was generalized to 
include the perturbation from the vector field in Ref. [5|. Up to the second order it reads 

C = N^Sdp + N^5W, + N^SW.SW, , (4) 

where iV^ = dN/dcj), N^ = dN/dWi and N^ = d^N/dWidWj and derivatives are with respect to the homogeneous 
values of fields. Also, the summation over repeated spatial indices is assumed. In this equation we also assumed 
that the scalar field perturbation is Gaussian, hence no second order terms in dcj), and that scalar and vector field 
perturbations are uncoupled. 

Using the SN formula we can easily find the spectrum of C at tree level i*^ 



Vc (k) = Pf (fc) 



i + ,c(k-N^)' 



(5) 



where fc = |k|. k = h/k and N^ = N^^ /Nw with Nw = |N^|. In this equation the isotropic part of the spectrum is 

7'r(fc)=iV|P^(A:)(f + e), (6) 



where Q 



-(^)^ 



and it quantifies the contribution of the vector field perturbation to the total curvature perturbation. For ^ < f the 
scalar field contribution dominates. 

The amplitude g^ of the quadrupole modulation in the spectrum in Eq. ([5]) is given by [5| 



^^V\\{k)-V+{k) e 



From the last equation we see that a mildly statistically anisotropic vector field perturbation, \p\ <C 1, can generate 
the total curvature perturbation, that is 5^ « p with i^ ^ I. Then the observational bound on g(^, discussed in 
the Introduction, gives p < 0.3. If, on the other hand, p violates this bound, the vector field can generate only a 
subdominant contribution to (^, i.e. ^ < f , in which case g,^ ~ ^P- 

Note, that the curvature perturbation power spectrum V(^ is proportional only to p but not q. Thus the possible 
parity violation of the vector field perturbation cannot be detected in the spectrum of C. We need to measure higher 
order correlators for that. 

In the present work we study the generation of ^ by the vector curvaton scenario, in which N^ and N^ can be 
found to be @ 



^"-^^-^ -d N-^^ICl^^^, (9) 



with Clw defined as 



In the above pw and p^ are energy densities of the vector field and radiation at the curvaton decay and 

Ow = pw/p, (11) 

where p — pw + Pr- If the curvaton decays while subdominant then Q.w — 3r2w/4 < 1. 



^ The quadrupole modulation of Vc^ can also be generated during weakly anisotropic inflationary expansion 13511 or inflationary models in 

non-commutative space-times [3|i]. 
* Note that this definition of ^ reduces to /3 in Ref. Q when V^ = V^. In this limit it is also equal to § in Refs. |19|| . 



The non-linearity parameters /nl for this scenario were calculated in Ref. f(| . In this reference two shapes of /nl 
were considered: the equilateral one f^, with ki k, k2 ~ k^,, and the squeezed one f^^-, with fci sa /c2 3> fc3. For flat 
perturbation spectra they are given by 



6 ,eql _ e 

r-'NL ~ /-, , ^^2 



{l+O 2^W 



i + i,' 



1 , 1 2 
P H P" Q 



Wf 



(12) 



fsqz 

'NL 



e 



{l + O 2^w 



1+pWl+ipqW^Jl 



Wf^ sin uj 



(13) 



In these expressions Wj^ is the projection vector of W = ~W /W onto the plane of the three wave- vectors ki, k2 and 
k3 with < W±_ < 1. The angle uj in the second equation is between Wj^ and the squeezed k vector. Here, we also 
calculate the third, so called flattened, shape where the length of one of the wave- vectors is twice as large as of the 
other two, e.g. ki k, k2 ~ ^k^ or ifci ss /c2 ~ ^3^ 



/■fit 
/nl 



e 



(l-fO 20 



w 



q' + [2p + p 



P 



cos (fiWi 



(14) 



where ip is the angle between W_l and the longest k vector. 

For the vector field contribution to C with g^ < 0.1 and p < 0(1) current observation constraints give flw > 10^^, 
with |/nl| ^ 100 [231 ■ The parity violation of the vector field perturbation (non-zero q) modulates the shape of /nl 
by suppressing its value in the equilateral configuration and enhancing it in the flattened one. It also introduces an 
imaginary term in the squeezed conflguration f^, which is real in position space by the virtue of the reality condition. 
From Eqs. P^ - plj) it is also clear that in general the vector field contribution to ( generates an angular modulation 
of /nl (terms proportional to W±). Large enough statistical anisotropy of the perturbation (parametrized by \p\ and 
\q\) can generate predominantly anisotropic /nl with a configuration dependent amplitude and form of the angular 
modulation. 



III. PARITY VIOLATING VECTOR FIELD 



A. Equations of Motion 



Let us consider the Lagrangian of a massive U (!) vector field 



(15) 



where Fny = dnA^ - di,Aa and 



1 pM'^PO' 



(16) 



with e^^'^P'^ being the totally antisymmetric tensor. The three functions /(t), h{t) and m? {t) are time dependent. 
Their variation is provided by other dynamical degrees of freedom in the theory. In this section we do not specify 
these degrees of freedom. Our aim is rather to find a scaling of these functions which give a flat perturbation spectrum 
for the vector field. 

The Lagrangian of the form in Eq. ([T5|) . with h — Q and non-zero mass term, was first proposed in Ref. [13] for 
the generation of the primordial curvature perturbation and then extensively studied in Refs. [19j. It was found that 
during inflation all three degrees of freedom of the massive vector field acquire a flat perturbation spectrum if the 
kinetic function / and the mass m evolves as 



/oc a 



-1±3 



and 



(17) 



where a is a scale factor. In this section we study perturbations of the parity violating vector field, 
evolution of functions / (t), h (i) and the mass m (t) which result in a flat spectrum. 



We find the 



The Euler-Lagrange equation with the Lagrangian in Eq. P^ gives the field equation for the vector field as 



The temporal and spatial components of this equation read as 



( nrn ] 

d,A, - d,d,Ao + ^^Ao = 



m^A" = 0. 



and 



A,. 



f 



H+Y\A, + ^A, - a-2 {d,d,A, - d.djA,) ~ a-^-e'^'d,Ai = 9, 



An 



"n '° 



While the integrability condition gives 

3iJ (d,d,Ao - d,A 
Combining these three equations we find 



(am) / m . ; o „ . 

^-^ [2-Ao + Ao-a-^d,A, 



A, 



H 



f 



A, - a-'d,djA, 



f 



-A 



a-^-f'^'d,Ai = ( ^ - 2- - 2if I d,Ao. 



f 



I 



(18) 



(19) 



(20) 



(21) 



(22) 



One can immediately notice that axial term changes only the equations of motion for the spatial components (and 
not the temporal one) of the vector field by adding the last term on the LHS. Also, as this term is proportional to the 
derivative of the vector field, the axial term does not have any effect on the homogeneous values of the vector field 
components. 

As noted in Refs. [l|, [131 the vector field A^, which enters the Lagrangian, is defined with respect to the comoving 
coordinates. While the physical vector field is (Aq, A^/a) for our choice of the metric. Thus, spatial components of 
the canonically normalized, physical vector field are given by 



W^ = V7— ■ (23) 

a 

To study perturbations 5Wi of Wi we go to the momentum space and write Fourier modes of 5Wi as 

SW,{\^,t) = e^{i^wx{k,t), (24) 

where summation over A = 'L', 'R' and '||' is assumed. Vectors e'", e^ and e" are three circular polarization vectors 
with k • e'" = k • e^ = k X ell = 0, k • ell = k and k x e'" = —ike^, k x e^ = ike^. Using this and Eq. (pi)) one easily 
notices that the axial term does not affect the longitudinal component of the perturbation. The equation of motion 
of this component is exactly the same as studied in Ref. ^19i] , where it is shown that a flat perturbation spectrum is 
achieved if Eq. (|17p holds, giving 



3H\- / H^ 



M 



In this equation M is the effective mass of the vector field 



PP 



27r 



,3±3 



/ 



(25) 



(26) 



where M = constant for f o^ a? . 

The transverse polarizations, however, are changed by the axial term. Modified equations of motion for the physical, 
canonically normalised vector field transverse modes are given by 



«^R,L + 3i?WR,L 



Af2± 



kh 

«7 



«^R,L = 0, 



(27) 



where wr,l {k,t) are functions of a wave-number k and we also used Eq. P7)) . The plus sign in the last term of this 
equation is for the right-handed polarization. Let us denote 



Q 



2 _ 



k 



a f 



Then Eq. ^7} can be written as 



w± + 3Hw± 



M^ ± Q2 



'w± = 0. 



(28) 



(29) 



If h (t) is positive during inflation, subscripts of w are understood as '+' = R and '— ' =^ L. If, instead, h (t) is negative, 
then '-I-' = L and '-' = R. 

It is difficult to find a general solution of Eq. (1^9)) . Thus we solve this equation in the regimes where each of the 
term proportional to w± is dominant and then match these solutions together. 

B. Initial conditions 

Initial conditions for each perturbation mode are set by assuming that deep within the horizon the k/a term 
in Eq. (|29p dominates. Effectively this means that such a mode describes a free quantum field and one can set 
Bunch-Davies initial conditions (see e.g. Ref. |37| ) 



ik/aH 



2k 



It is easy to see that with the above initial conditions modes w± and their derivatives are of the form 
where "H*-^-* are Hankel functions of the first kind. 



(30) 

(31) 
(32) 



C. Vector Field Perturbation Spectrum 

Solutions of Eq. (|29|) with the dominant M^ term can be found in Ref. ^9]. One could also consider the case 
Q ^ M. But cosmological scales span seven orders of magnitude in k. Because Q oc Vk, Q ^ M can apply only to a 
very narrow range of scales of interest, and we do not consider this case. In this section, instead, we study a situation 
when the dominant term is Q^, i.e. 

w± + 3Hw± ± Q^w± = 0. (33) 

For this we assumed the power law ansatz for h allowing us to write 

Q ex a". (34) 

A general solution of a second order differential equation has two constants. We determine these constants by requiring 
that the mode functions and their derivatives in Eqs. (|3ip and ([5^ match to the solutions of Eq. (155)) when the scale 
factor is a = flx, where 



Q (ax, k) . 



ax(fc) 
If Q is a power law function, Eq. (|33l) is a Bessel equation. For the w+ mode the solution of this equation is 

Q \ , ^+.. f Q 



w+ = a-3/2 



(-'I Jiy 



\c\H 



6*2 Yi, 



w+ = -a-^l'^Q 



L-i Ju+Sc 



Q 

\c\H 



c\H 

^2 ^i^+Sc 



Q 

\c\H 



(35) 

(36) 

(37) 



where J^ and Y^, are Bessel functions of the first and second kind respectively of order 



3 



(38) 



and Sc is a signature of c, i.e. Sc = ±1 if Q is an increasing or decreasing function respectively. By matching these 
solutions to Eqs. pTjl and ((32)) we find the constants C^ and C2 



+ _ 1 [T^T~Ll/2Yy — "^3/2^1/4 



1 2\ H J^Y^+ 



C^ 



1 /tt "^1/2 Ji/ — "^3/2^1/4 



(39) 
(40) 



where for brevity we used the following notation: Ti,^ = Hi, {k/a^H), J^ = J^ [k/ \c\ a^H), i/+ ^ ly + Sc and similarly 
for Y^. 

In the same way we can find the solutions for the w- mode 



w^ = a-3/2 



G]^ Xj/ 



Q 

\c\H 



G2 A^jy 



Q 



^ .«..,-3/2 



Q 



^l -Lu+Sc 



Q 

\c\H 



c\H 
— G2 r^v+Sc 



Q 

\c\H 



(41) 
(42) 



where T^, and ICi, are hyperbolic Bessel functions. The constants C^ and C2 are found again by matching the above 
solutions with Eqs. ^Q and ^^ 



„_ _ 1 rifScV. 1/2^1, + 'H'i/2^u+ 



(43) 
(44) 



where X^ = X^, {k/ \c\ a^H) and the same for K,^. 

The perturbation power spectrum for each mode is given by 



Vu 






(45) 



In the k/ucH — > limit we consider two possible values of the function Q: either Qc/H <C 1 or Qc/H ^ 1, where 
subscript 'e' denotes values at the end of inflation. In each case, the function Q can be decreasing or increasing. Thus 
for each value of Qc we have two possible values of fc/ox - larger and smaller than Qc - giving us four possibilities in 
total (see Fig.[T]). We study these four possibilities below and find the perturbation power spectrum in each case. 



The Case I 



This corresponds to the case when the inflationary Hubble parameter is represented by the uppermost dotted line 
in Fig. [2 i.e. k/a^H <C 1 and Qc/H <C 1. In this regime the Q term can be increasing as well as decreasing, i.e. 
c > or c < 0. We find the constants C^ and C^ by expanding Eqs. ([5^ and PO)) for small values of k/a^H 



IS 



Ct ^~^V{v + l) 



3 

iSr TT 



k 



2a^H^ \2\c\a^H 
k 



C+ = - 
2 3 r(zy) V 2axi?2 y2\c\a^H 



where 



l-[3+|c|(s,-l)]'^ 



k 



„+ — 



l-[3+|c|(Se + l)] '^ 



k 



a^H 



-l + Sc 



(46) 
(47) 

(48) 
(49) 




Figure 1. A log- log graph illustrating a possible evolution of the physical momentum kja and the function Q oc a'^ during 
inflation (for a fixed k). Horizontal lines represent four possible cases for the value of the Hubble parameter. The end of 
inflation at a = Be is depicted by a vertical dot-dashed line. The mode k exits the horizon when the blue (solid) line crosses 
the dotted line. The four possible cases are the following. In the Case I the function Q <^ H during inflation, no matter if it is 
increasing of decreasing. The Case II is when a function Q ^ H before the mode exits the horizon, but becomes smaller than 
H towards the end of inflation. This can only happen if Q is a decreasing function of time. The Case III depicts a situation, 
when Q <^ H before horizon exit, but Qe ^ H at the end of inflation. This only is possible if Q is an increasing function 
of time. Finally, in the Case IV Q ^ H for both increasing and decreasing Q during inflation. The black dots in the graph 
highlights the moment when k/a^H — Q. 



Note that, for an increasing Q, s^ ':^ §2 , and for the decreasing one s{^ ^ sj . In both cases, they are equal to 

s+ (c > 0) = s+ (c < 0) « -3 
Expanding the Bessel function in Eq. (J4T|) with Qd H <C 1 we find 



a^H 



-3/2 



ct 



c. 



T{v) ( Q 



T{v + l) \2\c\Hj ^ TT \2\c\H 

Inserting Eqs. pS)) . (IT7)) and taking into account Eq. (|50p . we finally obtain 



(50) 



(51) 



V2 



(52) 



Using Eq. pS)) we see that the spectrum of w+ is fiat 






(53) 



We apply the same method for calculating the power spectrum of the w- mode. In the limit of small k/a^H 
constants Cf and C2 in Eqs. (|43p and (|44| become 



where s^ ^-^d S2 are given by 



^i"=~^^(" + 'V^(^k^' ^^' 



^ 3r(i/) V 2a^H^ \2\c\a^H 



Si = Sc- [3+ \c\ (sc - l)Y 



S2 > 



a^H 



[3+|c|(sc + l)] "' 



a^H 



-l + Sc 



(54) 
(55) 

(56) 
(57) 
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One can see again that for c> 0, s^ ^ ^2 ^'^'^ '^^'^^ versa. Also for each of these cases 

-2 



Sj (c > 0) = S2 (c < 0) 
Next, expanding Eq. (PT|) for small Qc/ |c| i? we find 



'Sr3 



a^H 



-3/2 



Cf 



T{v + l) \2\c\H 



Qo 



C, 



T{u) ( Q 



2\c\H 



(58) 



(59) 



It is easy to find that, using Eq. (1581) and expressions for C^ and Cj in Eqs. (I54p and (|55p . the mode function w^ 
becomes 



V2 



Thus the power spectrum is also flat and equal to 



5(1) - 



' 27r 



H 



(60) 



(61) 



As we can see from Eqs. ([55)) and (pTj) . the power spectra for both transverse modes are identical at first order in 
Qe/H and are equal to the one of the light scalar field. This is not surprising as the equations of motion for w+ and 
w_ in Eq. ([55)1 reduces to the one of the light scalar field. The spectral tilt oiVw^, however, will differ from the scalar 
field case as the Q (t, k) is a function of time and momentum. 

To find the spectral tilt of Vw± we expand the Bessel functions to the second order in Qc/H . Let us consider Vw^ 
first. Then 






Q \ ^ 1 (_Q_ 

\c\h) T{l + v) \2\c\H 
Q \ T{v) f Q 



1 



\c\H 



11 \2\c\H 



1 - 



l + v \2\c\H 



Q 



l~v \2\c\H 



O 



O 



H 



v+4 



v+i 



(62) 
(63) 



The above expression for Y^ is valid only for v > 1. However, the Y,j term in Eq. (j36p is dominant only when this 
condition is satisfied. From the above, it is easy to infer that, to second order in Q^/ H ^ the power spectrum of w+ is 



^-- - ' I 



w+ 



1 



4c 



Qc 



3 + 2c \2\c\H 
Performing the same expansion for the hyperbolic Bessel functions we find 



X, 


' Q \ 


1 ( Q \ 


V 


[ll ' { "^ " 


2 

) 


+ 




T{l + v) \2\c\h) 


^ ' l + v \2\c\H^ 


JC 




r(^) / Q Y 

2 \2\c\h) 


1 


1 W "^ VI 


-I- ^-^ 


'fQY^^' 






' l~v \2\c\H ) 






\ 


Hj \ 





With this result, the power spectrum of w^ becomes 






1 



4c 



Qc 



3 + 2c \2\c\H 
The two expressions for Vw± above can be combined into one result 



(64) 

(65) 
(66) 

(67) 






^ 

'^'h 



(68) 
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where 

e= °'' A (69) 

'- c{3 + 2c)Hfo- ^^^^ 

As it is clear from the above, the spectral dependence of the power spectra of w^ and w- are the same but with 
opposite signs. The power spectrum of ^ in Eq. ([5]), on the other hand, is proportional only to the arithmetic average 
of Vw^ and Vw- hi which momentum dependence cancels out. Thus the vector field in this case may influence the 
spectral tilt of Vi^ only at even higher order in Qd H <^ 1, which probably makes it undetectably small. 

Spectral tilts of Vw± do not modify /nl significantly too. As can be seen from the definitions of the parameters p 
and q in Eq. ([3]) and the expressions for /nl in Eqs. (|T2|) - p4)) . the above spectral dependence infiuences /nl only 
through the parameter q. In this case, however, q (k) <C 1. 

2. The Case II 

Case II corresponds to when the Q term dominates the k/a one in Eq. ((29)) before horizon exit, but at the end of 
inflation Qe <C H. This is possible only if Q is a decreasing function of time, that is c < 0. As discussed in subsection 
IIIIB| however, to set initial Bunch-Davies conditions the k/a term in Eq. ([29]) must dominate initially. Therefore the 
Q term must decrease slower than k/a, putting the lower bound on c > — 1. 

In the limit k/a^H ^ 1, the constants C^ and C2 in Eqs. ([39| and pO|) become 



Ct^-,..Ct^^J-lj^^. (TO) 



where 



k f l\ 2v + l 

The equation for w+ in the limit Qc/H <C 1 is given in Eq. (|5ip . with the constants Ci and C2 determined by 
Eq. ([70]) . As V is positive by definition, the second term on the RHS of Eq. ([51]) dominates. Inserting it into Eq. (J45| . 
we find 



.,;.., ,,.-.=i^ay"m-7-)\ ,, 




One notices that with c=— 1/2, 1^ = 3 and the power spectrum of 'w+ is flat 




-3 . .2 



^SP = r(.;^) [^A • (73) 



Let's consider the w- mode next. In the limit k/a^H ^ 1 the constants Cf and C2 become 



and in the limit Qe/H ^ 1 the mode function w- is given in Eq. (|59p . Because Qc/H <C 1 and |C^| ^ I^TI the 
second term dominates and its contribution to the power spectrum is 

V2'} =V^f-e^''\ (76) 

where 

O / /, \ 2(c+l) / u \ 2T7TTT 

(77) 
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Here again one can easily notice that, with c = —1/2, the exponential amplification of the power spectrum ViuJ is 
scale independent 

em ^ 4 A. (78) 

3. The Case III 

Case III in Figure [T] corresponds to the situation when the Q term takes over the k/a term in Eq. (P^ after horizon 
exit, but Q is a growing function of time and becomes larger than H at the end of inflation. We find the power spectra 
of w+ and W- by the same method as in the previous two subsections. 

Since Case III corresponds to k/a^H ^ 1, constants C^ and C^ in this limit are given in Eqs. (|15|) and (|T7)l . While 
in the limit Qc/H ^ 1, Eq. ([55]) for w+ becomes 



-3/2 2|c|7f 



2i^ + 1 \ „+ . / Qo 2i/ + 1 



Ci+cos(^--^7r +C+sin 



-TT 



(79) 



For an increasing Q, i.e. c > 0, Cj^ 3> (^2 and the first term in the above equation dominates. Using this and 
Eqs. (gSl) we find 

where 

«'--^^-^- <-' 

As c > the power spectrum Vw^ in this case cannot be scale invariant. 

The integration constants of the W- mode for k/a^H <C 1 are given in Eqs. ([M]) and (|55p . Expanding Eq. ([1T|) in 
the limit Qc/H ^ 1 we find 



y 27r(5e ^ ^ 



(82) 



For increasing Q, the constants are C*]^ 3> Cj and the first term dominates. Using this we find 



T'i'^^^^i'^^Je^Q^/l^l^. (83) 

Because Qc/H 3> 1 is a function of a wave-number fc, we see that the power spectrum of u;_ features an exponential, 
scale dependent amplification factor. 

4. The Case IV 

In the Case IV, Q is always larger than the Hubble parameter during infiation. Such scenario corresponds to 
fc/ox ^ H and Qc ^ H. We have already calculated the integration constants and the equations for the mode 
functions in the previous subsections. The expression for w+ mode in the large Qc limit is given in Eq. ((79| . While 
the constants C^ and C^ in the large k/Ha^ limit can be found in Eq. ([70]) . With these equations it is easy to find 
the following solution for w+ 



„-3/2^+ /25^^...e<-)^ (84) 

y T^Qe 



13 

where the phase factor 0(m) is defined in Eq. ([5T|) . Substituting this into Eq. (P5)) . we find 

(85) 

As one can see this spectrum is very blue. 

In an analogous way, we can compute the power spectrum of the w- mode, which in the large Qo limit is given 
in Eq. ([5^ . The integration constants Cf and C^ were calculated in Eqs. ([7^ and ([75)) . Putting these expressions 
together we find 




a 

W- = 



-3/2 



2VQe 



g(Qo-fc/Qx)/|c|H _,_ g-(Qo-fc/a,)/|c|H 



(86) 



When Q is an increasing function of time, i.e. c > 0, then Qc ^ k/a^ and the first term in the above expression 
dominates. While, for c < 0, the second term dominates. With this solution, the spectrum of w^ becomes 

Thus, again in this case we find that the power spectrum of w- has an exponential fc-dependent amplification. 

5. Logarithmic h 

The last two cases - logarithmic h and constant Q - are not shown in Figure [TJ In the first, logarithmic h case the 
time derivative of h is constant during (quasi) de Sitter inflation: 



h = Hho ~ constant, 
where ho is the initial value of h. From the definition of the function Q in Eq. (|28D we notice that 

Q oc a-(l+")/^ (89) 



where a = — 1 ± 3 is the scaling of the kinetic function / ex a" . While from the definition of c in Eq. ([M| we can see 
that the logarithmic scaling of h corresponds to 

c = ±3/2, (90) 

where c = —3/2 is for a = 2. With negative c the function Q is decreasing. But in order to be able to impose the 
Bunch-Davies initial conditions Q has to be decreasing slower than k/a. Thus c = —3/2 is unviable. However, for 
a = —4, c = 3/2 and we can use the results of subsections IIII CT] IIII C 31 and IIII C 4l to deduce the power spectrum. 

6. Constant Q 

For the constant Q case, the equations of motion of both modes w+ and w^ in Eq. ([29| reduce to the familiar 
equation of a massive scalar field. The difference, however, is that the effective "mass" in this case is scale dependent, 
Q ex fc^ " and for the w_ mode this "mass" is tachyonic. For Q -^ H, the power spectrum is well known for this type 
of equation 

where 

- - /HI/- 

We see that in the limit Q -^ H, the spectrum is almost flat. It is also clear that both modes will have the same 
spectral tilt, but of the opposite sign: the spectrum of w^ is slightly blue-tilted and red-tilted for w_. However, as 
discussed in Subsection IIII C fj this does not affect the spectral tilt of P^ and makes a negligible contribution to /nl- 
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V^u + 


T^_ 


Case I 


{§r 


P.^ 


Case II 


mc\r-^^y{jiiyi^,y-\zr 




Case III 


(2|c|)^+^''£^(^^) ' (f)^cos^e<-) 


1 Wj^ ^ c 


Case IV 


\Hf^) [a^Hj {27V ) 


T3 1 2(Q<,-fc/ax)/cH 
/ w+ 2 '^ 


h cclna 


The same as cases I, III and IV with c — 3/2 


Q = constant 


V 4p2/ w£r\2/ fe \3- 


2,.± 


1 w± ^1 IZ-iJU^rj {2a,H ) 





Table I. Summary of subsection HITCl In this table e*"' is defined in Eq. (|73, 9'"'' is defined Eq. dHJ and i/± in Eq. ((92)l . 



7. Summary of the Subsection 



We summarize the results of subsection HIT CI in Table H] There are two parameter spaces for producing a flat 
perturbation spectrum for both of transverse modes. First, this can be realised for any value of c if the Q term in 
Eq. (P^ dominates the k/a term after modes exit the horizon and if Q stays smaller than H until the end of inflatiorO- 
Such a vector field can generate statistical anisotropy in the curvature perturbation C, which is consistent with the 
observational bounds. This can be realized, for example, using the vector curvaton scenario [l|. In this case the 
results of the vector curvaton scenario discussed in Ref. [1^ are directly applicable But as both transverse modes 
acquire the same perturbation amplitude, the axial term in the Lagrangian will not have a detectable signature in the 
curvature perturbation in this parameter space. 

However, in this paper we are also interested in a possibility of producing C, with parity violating statistics. As 
shown in Sec. IIII(T2] this can be realized if the function Q scales as c = —1/2 and if it dominates k/a term in Eq. (|^5)) 
before horizon exit. Note, however, that Q ex fc^'^. Therefore the above discussion is valid only for some limited 
range of k values. Particularly, Vw+ — V^^ — {H/2-k) if k/a^H <C 1 and Q^k) /H ^ 1. And parity violating 
perturbations are realized if the former bound is reversed and c = —1/2. We study the latter case in more detail in 
SecHVl 



IV. STATISTICALLY ANISOTROPIC, PARITY VIOLATING CURVATURE PERTURBATION 

A. The Spectrum 



In Subsection IIII CTj we found that the vector field acquires a scale invariant perturbation spectrum if the Q term in 
Eq. (|29p dominates the momentum term k/a after the mode k exits the horizon and Qc ^ H at the end of inflation. 
In this case both modes, w^ and ui_, acquire perturbation spectrum equal to (H/2t:) . While in such a set-up the 
scale invariant curvature perturbation can also be generated, it does not give any signature of parity violation. 

A more interesting case is studied in Subsection IIIIC 21 In this case, we find that the perturbation spectrum of the 
vector field is also scale invariant if the Q term dominates the momentum one in Eq. (j^^ before the mode exits the 
horizon and if Q is decaying during inflation as Q ex a~^/^ with Qc ^ H &t the end of inflation. From the definition 
of Q in Eq. (P5)l we see that this scaling of Q corresponds to 



h 



constant. 



(93) 



^ We assumed here M <^ Q. 

® The only difference are the additional constraints k/a^H <^ 1 and Qe/H <C 1. 
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With the power law evolution of / ex a" the above condition implies that h also evolves according to the same power 
law /i oc a" 13, where a = — 1 ± 3. In view of this, we can rewrite the Lagrangian in Eq. (jlSp as 



^—4 



\f (f^^uF^^'^ + ^F^.F^-') + ^m^A^Af^, (94) 



where 



d = — = constant. (95) 

Then a vector field with the Lagrangian in Eq. ([M)) will acquire a flat perturbation spectrum if the Q^ = |a^| Hk/a 
term dominates the (k/a) term in Eq. ([29| before the mode exits the horizon. In the case when Q ex a^^/^ this 
happens at a scale factor a^ given bjo 

i- h 

— - = \ad\ . (96) 



Ha^ Hf 

From this we see that the value of k/a^ does not depend neither on time nor on wave-number k. 

To calculate the perturbation spectrum of the vector field we used Bunch-Davies initial conditions in Subsec- 
tion IIIIBI These initial conditions are valid in the limit where the curvature of space-time can be neglected and 
quantum fields are effectively described by the quantum theory of free fields. The first requirement is fulfilled for 
modes deep within the horizon, while the standard quantum field theory can be applied for modes which are not too 
close to the Planck scale, and the mode can be considered to be of an effectively free field if the k/a term in Eq. ((29|) 
dominates. Using Eq. (|96|) we find that these requirements constrain the value of •& as 

1 < \a^\ « ^. (97) 

Because the parameter a is of order one, we see that the scale invariant vector field perturbation spectrum is 
achieved if the parity violating constant •& is much larger than unity. Then the spectra of modes w- and w^ from 
Eqs. ([73]) and ([76|) can be written as 



4, ..-3/H 



T'.^^-Mr^^j'; (98) 

Note, that due to Eq. (|57|) the spectrum of w+ is suppressed with respect to the standard light scalar field result. 

For the expressions in Eqs. (|M)) and (|M)) to be valid, the Q term must also be much smaller than H by the end of 
inflation, Q^ ^ H. It follows from Eq. ([M|) that this implies 

la^l « ^, (100) 

which must hold at least for cosmological scales. Assuming the observable inflation to last at least 50 e-folds, this 
bound is much weaker than the one in Eq. (1971) for any realistic model of inflation. 

If for the largest k modes this bound is violated, we are back to the Case IV in Figure [1] which was discussed in 
Subsection IIII C 4l As one can see in Eqs. ([55)) and ([57)) . the spectra of both modes are very blue Vw± oc A:^/^. As 
the Vw- is also exponentially enhanced by exp(4|ai?|) one may worry about the overproduction of the primordial 
black holes in such a scenario. The summary of various cosmological constraints on the abundance of the primordial 
black holes can be found in Ref. [3a|- On practically all scales the bound correspond to V( < 10^^ for the Gaussian 
perturbation. If the perturbation is non-Gaussian then the bound becomes 

-3 



n (^Pcak) < 10--^ or 1, (101) 



^ Remember that we assume (quasi) de Sitter inflation, i.e. H ^ constant. 
® The moment fc/ax is represented by black dots in Figure [T] 
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where the lower bound is for the positive non-Gaussianity and the upper bound is for the negative non-Gaussianity 
|39| . The contribution of the vector field to the spectrum of ( is N^Vw_, where N^;- is defined in Eq. ^. From 
GOBE normalization this contribution has to be < 10~^ for cosmological scales. From Eqs. ([98| and ([99|) Vw- ^ "Pw^ 
and thus the spectrum of C is 

N^-\a^r'(^ye'\'^'\<lO-', (102) 

for cosmological scales. 

The spectrum of modes which violate the bound in Eq. (|100p is given in Eq. ([57]) . The blue spectrum peaks at the 
largest k value, which corresponds to the horizon size at the end of inflation. Thus setting fepoak/o-e-ff ~ 1 and using 
Eq. l[87|) we can write the bound in Eq. (jlOip as 

iiV^ \a^\-'/^(^^y e^l-^l < 10-3 or 1, (103) 

where we also used the fact that fc/cx ^ Qe for c < 0. Using Eq. p02p this becomes 

|a^| < 10^ or 10"*. (104) 

As shown in Subsection (jIII A[) . the equation of motion for the longitudinal mode of the vector field perturbation 
is not affected by the axial term. Thus the results of Ref. [l^J for this mode can also be used directly in our case. It 
was found that the spectrum of the longitudinal mode is 

3HV / H^ ^ 



^'HjT S). (10=) 



where M is an effective mass of the vector field defined in Eq. ([261) and the index 'e' indicates that it is evaluated at 
the end of infiation. In these works it was also found that the spectrum of transverse modes are unaffected by the 
mass term if the field is light. 

The predominantly parity violating perturbation will be generated if Vw_ > Vh. Gomparing Eqs. (|99p and (|105p 
we find the condition for this to be 

Mr^^|"''|>(^)'. (106) 

Taking into account the bound in Eq. (|104p one sees that the above condition can be easily satisfied even for extremely 
small values of M^ . 

Results for the power spectrum given in Eqs. (|98l) and (j99l) are valid if the vector field is light during inflation. If 
this is not the case, and the field becomes heavy before the end of inflation, it starts oscillating and the amplitude of 
the spectrum decreases by i {3H/M) [l3- As can be seen from the definition of the effective mass M in Eq. (P5)l it 
can grow only if a = — 4. In view of this we can rewrite Eqs. ([M)) and (|99p as 

V^^^ — m-'f-X-mmll,—]'; (107) 

+ 167r ' ' \2ttJ 2 \ ' AfoJ 

B. The Vector Curvaton 

So far we have discussed the perturbation of the vector field. One way a vector field can generate or contribute to 
the primordial curvature perturbation ^ is via the vector curvaton scenario [l|. In this scenario it is assumed that 
the vector field is light during inflation, at least while the cosmological scales exit the horizon. To avoid excessive 
anisotropic expansion of the Universe, the energy density of the light vector field has to be subdominant. After 
infiation, when the vector field becomes heavy, it oscillates with a very high frequency and behaves as preasureless. 
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isotropic matter. The energy density of such matter decays slower than the radiation, and thus the vector field can 
dominate and generate the total of C or nearly dominate and generate a contribution to ^ 

C={i-nw)Ce.d + nw<:w, (109) 

where flw is defined in Eq. (|10|) and (^w is the vector field contribution to the curvature perturbation (. If the vector 
field perturbation spectrum is statistically anisotropic, so is the spectrum of ^ (SljlSl. 

In Ref. [l9| a vector curvaton scenario was studied with the Lagrangian in Eq. ([94| and ?? = 0. Bounds derived in 
these references also apply to the current scenario as the axial term with •& ^ does not contribute to the homogeneous 
values of fields and, therefore, to the homogeneous value of the energy density. This is because for the homogenised 
vector field, Fij — 0, where i,j denote spatial components. The axial term is proportional to e'"^'"^F;^,yFpo-- This 
means that, it cannot include any term featuring simultaneously two factors of the form Foi = —FiQ — A, which is 
the only non-zero component of the field strength tensor. Thus to find the constraint on the energy scale of inflation 
we can use the expression for the Hubble parameter found in Ref. ^191] 

H ,„ . r MA-'" . r ml^/^Vmax^Vl^WV^' mm 

— ^l^^CH^mmjl;— I mm|l;-| (^ j . (110) 

In this equation H is the inflationary Hubble parameter, fi-yv is defined in Eq. (jlip . T and Tw are the inflaton and 
vector field decay rates and i?dom is the value of the Hubble parameter after inflation, when the oscillating curvaton 
dominates over the radiation bath, if it does not decay earlier. 

With the axial term the dominant contribution to the vector field perturbation comes from the W- term, thus 
Vw^ > 'P«,+ , V\\ and V+ w jPw_ . Using Eq. ^ we find 



5C 



1 N^V^_ 

2 7^'^° 



(111) 



Also, using the fact that the anisotropic contribution to the power spectrum of C, must be subdominant [10] (cf. Eq. ([T]) ) , 
to the first order we can write 



/-2 -piso 



^ (112) 



5W^ ~ V^_ ' 

where we also used the fact that the dominant contribution to the vector field perturbation is from the w_ mode. 
Combining the above two equations we find 

N^SW"" ~ -2g^C- (113) 

Comparing the 5N formula in Eq. @ and the equation for C, in the curvaton scenario in Eq. (|109|) we can also write 

Nw5W ^ClwCw- (114) 

Inserting this into Eq. (|113p we get 

C-^^, (115) 



which can be used in Eq. (|110p to find the lower bound of the inflationary Hubble parameter. In Ref. ^3 it was 
shown that the lowest decay rate of the vector field is through the gravitational decay, which gives max {F^y ; ffdom} ^ 
M|?/m|j. Using this and Eq. (|115p we find the bound on the inflationary Hubble parameter H 

^>('-^V'7^V'7^V'''min/l-^r'''minIl>V (1161 

As we can see from this inequality the bound is maximised if the Universe undergoes prompt reheating, that is 
F ^- if*. Also note, that the bound on 7J* is maximised for the smallest value of the effective mass of the vector field 
Mq. Requiring that the vector field decays before the Big Bang Nucleosynthesis (BBN) gives Mc ^ 10'* GeV ^191] . 
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Inserting these limits into (|116|) and using the fact that i^w < 1, the bound on H^, and inflationary energy scale V^ 
becomes 



1/4 



H > V^IO^ GeV ^ V^/j'^ > (-.gc)^/'' 10^^ GeV. 
Also from the results in Ref. [19] the constraint on the mass of the vector field is 

10 TeV < A/e < 10^ H, 



(117) 



(118) 



where the bound in Eq. (|117p is saturated with the lower bound in Eq. (|118p . 

From the above two constraints it is clear that there is an ample parameter space for this scenario to be realised. 



C. The Parity Violating Non-Gaussianity 

As we saw in Section |ll] to calculate /nl it is convenient to use parameters p and q defined in Eq. ([S]). Because Vu 
is exponentially larger than both Vw+ and Vu , p and q are equal to 



p « —1 and \q\ w 1. 



(119) 



The sign of q is determined by the sign of a-d. Following the definitions of w+ and w_ after Eq. (P^)) we find that 
q w — 1 if a-d > and g sa +1 otherwise. 

The value of \p\ « 1 violates observational bounds on the anisotropy in the spectrum. Thus, as discussed in 
Section pl| . the vector field contribution to C, must be subdominant, i.e. ^ < 1, which gives 



5C 



-^ and Clw ~ -7^w- 



(120) 



Putting the above values of p, q, ^ and Clw in Eqs. p^ - P^ we find 



6 „eql 

5^' 



3^ 



W \ 4 

2 



-,2 



?/nT = 2^ (l - ^1 - * sgn{ad)W^.JT~Wlsh 



la = ^1^ {I- cos' ^Wl) 

O llw 



(121) 
(122) 
(123) 



where sgn(ai?) is the sign of a-d. As Vw^ ^ ^|| , the above result is valid for both: the light vector field and the one 
which becomes heavy at the end of infiation. 

The shape of /nl given in Eqs. (|121l) - (ll23p provides a smoking-gun signature for this model. It is easy to see that 
the maximum values of /nl in different configurations are related as 



1 roql 

o^/nl 



Re 



rsqz 
JNL 



5 fflt 
^/nl 



(124) 



while only f^ has a non-vanishing minimum value of 3/4. 



V. CONCRETE EXAMPLES 



In this section we consider some specific models motivated by particle physics, which can provide a flat vector field 
perturbation spectrum with the use of the axial term ex FF. But before going into the models it is important to 
point out again that the axial term does not affect the dynamics of the homogenised fields. Thus, we will ignore this 
term when discussing the dynamics of the homogenised scalar and vector fields. 
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A. String inspired model 

In string theory the gauge kinetic function and the couphng of the axial term are of the following form: / = RcJ^ 
and h — ImJ-", where J- is some complex holomorphic function of the moduli fields. Thus, the gauge field content of 
the model is 

C = -i(Rc^)F^,F^'^ + im^A^A'' + C (ImJ■)F^,F^^ (125) 

The simplest form for J^ is J^ ex e^^^ . for large values of the modulus T in Planck mass units {h = constant). Writing 
T in terms of the real fields (f) and a a,s T = (f> + ia and reinstating the Planck mass we find 

/ = e~'"^/™fi cos(6cr/mpi) 1 

I ^ - =^9 = Ctan(6cr/mpi), (126) 

h = Ce-'"^/'"Pi sin(6cr/mpi) J 

where x9 was defined in Eq. (j95p . If the axion a remains frozen during inflation with a ~ constant then we have 
"d = constant, i.e. / oc /i oc a" and c = — ^ This is a very realistic possibility because the axion mass is protected 
by the approximate U(l) symmetry and can be very small, much smaller than H during inflation. Note that we 
have not clarified whether (j) and a are canonically normalised fields. Indeed, if the Kiihler potential has a non-trivial 
dependence on these fields they will not be canonically normalised. However, this does not change the result above, 
that is c = — 2 • 

In the opposite limit of small T the dependence of T on the modulus is logarithmic. For toroidal compactifications 
the Kahler potential is of the form K = —3 ln(T + T) in Planck units, where T is a complex structure modulus here. 
This means that K = — 3mpjln(20/mpi). Hence, the kinetic term of the (j> field is non-canonical and it reads 



V 

In view of this we can define the canonically normalised scalar field $ as ln{(j>/mpi) = A=^/mpi. Thus, if we ignore the 

axion, the gauge kinetic function is / ex ReJ^ ex (ln0)" ex $", i.e. it has a power-law dependence on the canonically 
normalised field <&. 

In summary, provided we ignore the axion (presumed frozen), string theory suggests that h/f = constant. The 
functional dependence of / and h on the varying modulus (which could be the inflaton) depends on the (unknown) 
compactification scheme but both exponential and power-law dependence is reasonable. 

The condition k/a^ ^ H is ensured by the lower bound in Eq. ([W)) . i.e. -d ^ 1/4 with a = —4. From Eq. p26p . 
assuming C '^ 1, we see that this condition translates into the lower bound for a 

-^>10"V&- (127) 

mpi 

Thus, assuming 6^1, this condition is satisfied provided the value of the (frozen) axion is comparable to the 
Planck scale, but this can be relaxed if C ^ 1 (see for example Ref. [29]). The upper bound on ■&, which limits the 
overproduction of primordial black holes and is given in Eq. (|104p , translates into 

Ctanf -^ ) < lO^orlO^. (128) 

\mpij 

As was mentioned before, this also ensures that Qc -^ H for cosmological scales, as the latter is a weaker constraint. 
If the two bounds in Eqs. (|127p and (|128p are satisfied, it is possible in this set up to generate parity violating ( with 
a flat power spectrum, as discussed in Section HVl 

B. The orthogonal axion 



This model is of the form 



£ = D^^D^'S?)* - y($) - i/^^.F^'^ + ce^0F,,,Ff''' 

= ^d^<l>d^cf, + ^d^adf^a - y(0) - V{a) - 

--JF^^F^- -f l-e^cf'^A^A^ + ce2^F^,i?^^ (129) 

4 2 00 
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where D^ = 9^ + ieAf^ and / is modulated by some inflaton field but </> or cr are not it. Comparing the above 
Lagrangian with Eq. p3|) we see that 



h^~4ce^ — 



(130) 



in this model. To have some time-dependence for the axion a = 6(J)q we need its mass to be comparable to the Hubble 
parameter. Thus, we consider modular inflation with H ^ TeV and we assume that a is an axion field, orthogonal 
to the QCD axion in supersymmetric realisations of the Peccei-Quinn symmetry, which employ a non-renormalisable 
superpotential for the Peccei-Quinn fields |40l l4l|. This construction also solves the ^-problem of supersymmetry. 
The superpotential is of the form 



W 



K $"+3 



n + 3 m 



PI 



where $ — {(j)/^/2) exp{i9 / ^/2) , with 9 = a /<J)q with 00 being the Peccei-Quinn breaking scale given by 



= 2 



„ \ l/(n+l) 



V2k 



(131) 



(132) 



where m^ ~ TeV is the tachyonic soft mass of the radial field 0, which breaks the Peccei-Quinn symmetry. The above 
is obtained by minimising 



1 2 A2(n+2) 



(133) 



where Vq is some density scale. The orthogonal axion potential is 

y(a) = (0m,)2[l~cos(a/(/.)]. (134) 

The above potential is due to the soft A-term in the scalar potentially which gives the following value for the mass of 
the orthogonal axion: 



mi — kA 



(0/v^)"+i 



m 



(135) 



PI 



with A ^ TeV being the coefficient of the A-term. If the radial field assumes its vacuum value (p = (j>o (i.e. it is not 
rolling down the radial direction) then the axion mass becomes 



ll = 2("-2)/2^ 



(136) 



We assume that, during inflation, the axion is very close to the origin and rolls down near a local minimum 
approaching the origin. Then we can write 



dVia) 



da 
Using the above, the axion's equation of motion is 



™^'/'osin(cr/(?!)o) ~m^cr. 



a -I- "ilia + m„a ~ . 



whose growing mode solution is 



a ex exp ■ 



1-Wl 



2 ?7V 

sir 



Ht y oc 



„-ih^Hf¥r], 



(137) 



(138) 



(139) 



^ 5Va— term = A{W + W*) for a monomial superpotential. 



Thus, if we denote 



from Eq. P^ we find 
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(140) 



1-Wl 



2 7TV 

3 H 




2 m^ 
3~h" 



>0, 



(141) 



where we also considered Eqs. ([28]) and ^^ which suggest 2c == 7 + 3. If mg. <C i? then 7 « ^i-^)^ < 1 ^ c ~ 3/2, 
i.e. (T slow-rolls (and h oclna). In fact, a is practically frozen, in which case there is no parity violation. For this 
reason we do not consider that the radial field is still rolling with </> <C i^Oj neither do we assume that the tachyonic 



mass of the radial field is suppressed by some supergravity correction so that (m^)eff = rrij 



cH"^ < ml 



1 TeV, as 



H ^ 1 TeV, according to 



in Ref. [41]. Both these possibilities would result in ma <C -ff ~^ 1 TeV. Therefore, we should assume that the Peccei- 
Quinn symmetry is fully broken with (j) — (J)q and to^ ^ 1 TeV, such that rria 
Eq. (fnSD . 

The above imply that the mass of the physical vector field is 



v^ 



M 



— = = — = oc a . 

Vf Vf 



(142) 



Thus, the longitudinal component of the vector field exists and, if it undergoes particle production, it will not obtain 
a scale-invariant spectrum. We need to check whether the longitudinal component spoils the model. 
From Eqs. P^ and (fTi^ we find 



M 
ll 



m3/2 



n/{n+l) 



-2^ ^ e X 10i5"/("+i)e- 



■2N 



(143) 



where 7713/2 ^ 1 TeV stands for the weak scale (gravitino masa h and we considered that, at the end of inflation, 
/ — ^ 1. Firstly, we need to verify that the field is light when the cosmological scales exit the horizon. The vector field 
becomes heavy (and begins oscillating) when M ^ H. The earliest time for this to happen can be found by taking 
e = 1 and n -^ 00 in the above. We obtain N"^^^ ~ -^ In 10 ~ 17, which corresponds to much later times than the 
exit of the cosmological scalesFH 

Secondly, we need to ascertain that the vector field, after it begins oscillating, survives until the end of inflation 
at least, so that it can have some hope to affect the curvature perturbation in the Universe. The decay rate of the 
vector field is Tw — ^M oc a . Requiring that Tw ^ H at the end of inflation produces the constraint 



< 3 X io-W("+i) 



0.01 for n = 1 



(144) 



10- 



for n 



CX) . 



This is a tight constraint for the gauge coupling and excludes large values of n. Using the above, we can estimate 
how close to the end of inflation the oscillations of the vector field begin. Indeed, it is straightforward to find 



N„, 



< 



lOn 



In 10. 



6 for n = 1 



1 1 for n — > 00 . 



(145) 



Thus, in all cases the field becomes heavy and oscillates a few e-folds before the end of inflation, but it is always light 
when the cosmological scales exit the horizon. 

Now, let us consider the particle production process for the longitudinal component. In Ref. [19| it was found 



T'li 



167r 



sin2(7ri>)[r(l-z>)] 



hV f hV f k 



27r 



Mj \2aH 



b-2v 



(146) 



for gravity mediated supersymmetry breaking. 



N* 



1/4 



50 for prompt reheating with V-^^f ~ ^mpim3/2 
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where 

i'=^^/9 + 2{a + 1)(2 - a + 2/J) + (2 - a + 2^)^ , (147) 

with f (X a"^ and m^ ex a*^. Using that, in this case, a — —4 and /3 = we find i) = ^, which gives 

" \2ttJ \MJ \aH ) V " 2-11 aM ^ ' 

As we have discussed, at the end of inflation, M ^ H. Thus, since we are considering superhorizon scales with 
k <^ aH, we see that, at the end of inflation Vw ^ V^ = {H/2tt)'^. Hence, even though Vw ex. fc^, the scale invariance 
of the curvature perturbation is not spoilt because the contribution from the longitudinal component of the vector 
field is negligible compared with the transverse components, which are (approximately) scale- invariant. 

In this model flat perturbation spectra for transverse modes will be produced if fc/ox ^ H and Qc <^ H (c.f. 

subsection IIII C l|) . One can see from the definition of Q in Eq. (pS)) that {k/a^H) — ( h / f) = 'y{h/f)\^. Now, 

from Eq. pUjl with m„ ^ H ^ 1 TeV, we get 7 = 0(1). We also find h^ = ce^{cr^/(j)o) < 1, since c = 0(1), e < 1 
and < (Tx ^ T^<Po- Finally, because /x — e*^" 3> 1, the physical momentum at k/a^ — Qx is k/a^x. ^ H as desired 
(iVx here denotes the remaining infiationary e- folds when a — ax)- Similarly, Q^ is also smaller than H. Since at the 

end of inflation, on superhorizon scales k ^ {(^H)^ and CTo ^ (t>o '^e find [Qd H) = (^) ( f^ ) ^1- Thus, we see 
that we do obtain scale-invariant spectra for the perturbations, which means that this model can be used to generate 
statistical anisotropy in C,. However, as discussed in the end of Sec. IIIIC H in this case, parity violating signatures 
are probably undetectable both in the power spectrum and in the bispectrum, whatever the value of c. 

VI. CONCLUSIONS 

In conclusion, we have demonstrated that it is, in general, possible to employ an axial coupling for a vector field 
in order to generate scale invariant statistical anisotropy. There are two possibilities for this. In the first possibility, 
the only requirement is that, after horizon exit, the contribution of the axial term dominates the momentum term in 
the equation of motion for the transverse vector field components (the longitudinal one, if it exists, is not affected by 
the presence of the axial term) but still remains smaller than the Hubble scale until the end of infiation. In this case, 
both transverse components obtain a scale invariant, superhorizon spectrum of perturbations of magnitude i?/27r. If 
the field is massless then there is no longitudinal component (it is decoupled from the theory) so the anisotropy in 
the particle production is 100%, which means that the vector field has to contribute subdominantly to the curvature 
perturbation C,, as discussed in Ref. [5]. If there is non-zero mass, the longitudinal component is physical and has to be 
taken into account, in the way described in Ref. [131 • The vector field can play the role of vector curvaton and produce 
statistical anisotropy in the spectrum and bispectrum of C, as in Ref. [13- It can also contribute to C via another 
mechanism, e.g. the end of inflation mechanism [J,l5|. We have presented one example based on particle physics, where 
the axial coupling involves the so-called ort hog onal axion, which is orthogonal to the QCD axion in supersymmetric 
realisations of the Peccei-Quinn symmetry [40, |41|. In this case, the axial contribution to the equations of motion is 
growing during inflation, while there is a longitudinal component which obtains a scale-dependent spectrum, that is 
subdominant to the scale invariant spectrum of the transverse components. However, this possibility does not produce 
any parity violating signatures in the bispectrum. 

In contrast, the second possibility can indeed generate a parity violating signature. This possibility corresponds to 
a decreasing contribution of the axial coupling to the equation of motion of the transverse vector field mode functions. 
This contribution has to take over the momentum term before horizon exit. The scenario can be realised only when the 
couplings / and h of the kinetic and the axial terms respectively are proportional to each other, with proportionality 
constant d = h/ f . Then the power spectra of the transverse components are 

7?^ =l|a^|-3(iy/27r)2 and V^ ^-\ad\-'^e^\''-^\{H/2nf, 
TT n 

i.e. Vw^ — ^e^^"'^^'Pw+, where a — —1 ±3. As we have discussed, this possibility can be naturally realised in string 
theory, where / and h are determined by moduli fields. The parity invariant signature will affect the non-Gaussianity 
by modulating the angular dependence of /nl- If the Planck satellite does detect statistical anisotropy and anisotropic 
non-Gaussianity comparison between the different configurations (e.g. equilateral, squeezed and fiattened) may well 
reveal the existence of a parity violating signal and provide evidence of an axial coupling in the vector field which will 
be needed to explain the statistical anisotropy. 
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